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Abstract 

We demonstrate that lattice QCD calculations can be made 10 3 - 
10 6 times faster by using very coarse lattices. To obtain accurate 
results, we replace the standard lattice actions by perturbatively- 
improved actions with tadpole-improved correction terms that remove 
the leading errors due to the lattice. To illustrate the power 
of this approach, we calculate the static-quark potential, and the 
charmonium spectrum and wavefunctions using a desktop computer. 
We obtain accurate results that are independent of the lattice spacing 
and agree well with experiment. 







1 Introduction 



For 20 years, quantum chromodynamics (QCD) has been the generally 
accepted theory of strong interactions in particle physics. However, physicists 
are only just beginning to show that it explains basic features of low-energy 
strong-interaction physics, such as the spectrum and structure of the hadrons. 
The most significant progress in the study of low-energy (nonperturbative) 
QCD has come from numerical simulations of a lattice approximation to the 
theory. However this approach has been severely limited by the rapid rise 
in computational difficulty as one approaches the physically relevant limit of 
small lattice spacing and large lattice volume. In particular the cost of lattice 
simulations typically grows as a~ 6 when the lattice spacing a is reduced. 
With standard lattice techniques, it is widely felt that lattice spacings of 
0.05-0.1 fm or less are necessary for reasonable accuracy, making simulations 
impossible on all but the largest computers. In this paper, we suggest 
modifications of lattice QCD that permit accurate simulations at lattice 
spacings as large as 0.4 fm, increasing the speed of simulations by a factor 
of 10 3 -10 6 . We illustrate our approach with state-of-the-art calculations of 
the static-quark potential and the charmonium spectrum performed on a 
high-performance personal computer M. 

The central approximation in lattice QCD is replacing continuous 
spacetime by a discrete lattice. The QCD action is discretized by replacing 
space-time integrals with sums and derivatives with differences. Then 
the path integral defining the field theory can be evaluated numerically 
using Monte Carlo techniques. Standard discretizations of the QCD action 
have errors of 0(a 2 ) that are large when the lattice spacing is 0.4 fm. 
However improved discretizations can be designed in which finite-a errors are 
systematically removed by introducing new (nonrenormalizable) interactions 
into the lattice action. This does not add new parameters to the theory, since 
the coefficients of the new interactions are determined, from first principles, 
by demanding that the discretized action reproduces continuum physics to 
a given accuracy. Since the new interaction terms correct for deficiencies 
in the short-distance behavior of the lattice theory, their coefficients can be 
computed using perturbation theory in asymptotically free theories such as 
QCD, provided the lattice spacing is small enough that perturbation theory 
is applicable at distances of order a and smaller. 

The use of perturbatively improved actions for lattice QCD was suggested 
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long ago [[J . However early tests showed little benefit from the improvements. 
Furthermore a variety of studies seemed to indicate that perturbative 
dynamics was only relevant at distances significantly smaller than 0.1 fm, 
suggesting that lattice spacings would have to be smaller as well. Recent 
work || on lattice perturbation theory, however, has completely changed 
this conclusion. If conventional lattice perturbation theory (in terms of 
the bare coupling) is replaced by an expansion in terms of a renormalized 
coupling, perturbation theory becomes useful even at distances as large 
as 0.5 fm. If in addition lattice operators are "tadpole improved" the 
convergence of the perturbative expansions needed to define improved actions 
is greatly enhanced. Indeed, as we show here, tree-level calculations of the 
improvements suffice in most cases provided operators are tadpole- improved. 
Without tadpole improvement, the correction terms are systematically 
underestimated, sometimes by factors of two or three. 

Tadpole improvement is a technique for summing to all orders the large 
perturbative contributions that arise from tadpole diagrams peculiar to 
lattice QCD. At tree level the improvement is trivial to implement: one 
replaces each lattice QCD link operator in the action by = U^/uq, 
where «o is a scalar mean value of the link, defined to be the fourth root of the 
expectation value of the four-link plaquette in Monte-Carlo simulations Q. 

In several earlier papers, we have advocated the use of perturbatively 
improved actions with tadpole improvement |], |3|, £|, |j| . Such actions have 
already proven very successful in simulations of heavy-quark mesons like 
the T. In the nonrelativistic NRQCD quark action ||, both relativistic effects 
and finite-a corrections are introduced through nonrenormalizable corrections 
to the basic action. The detailed simulation results presented in || agree well 
with experiment, and many depend crucially on these corrections. Tadpole 
improvement was essential to this success; test simulations without tadpole 
improvement underestimated relativistic effects by as much as a factor of 
two. Similar results have been obtained with improved versions of the Wilson 
quark action when applied to heavy-quark mesons [|7| . 

In this paper we present new evidence that perturbatively-improved 
actions, once they are tadpole improved, work well for gluons as well, even 
at spacings as large as a = 0.4 fm. Using improved actions for the quarks 
and gluons on a lattice with a = 0.4 fm, we obtain a static potential 
that is rotationally invariant to within a few percent, the spin-averaged 
charmonium spectrum accurate to within 30-40 MeV, and rotationally 
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invariant charmonium wavefunctions. Our results show lattice-spacing 
independence (scaling) to within 3%. 



where 



2 The Improved Action 

The standard Wilson action for gluons has finite-a errors of order a 2 . On 
coarse lattices these lattice artifacts lead to severe (up to 40%) deviations 
from rotational invariance of the static quark potential. This can be seen 
clearly in Fig. |TJa (the potential computed on a lattice with a « 0.4 fm), where 
the points for r/a — \pl and v^3 lie far off the line defined by r/a — 1, 2, 3. 
These artifacts arise because the plaquette operator, from which the Wilson 
action is constructed, contains 0(a 2 ) terms beyond the desired gluon kinetic 
term when it is expanded in powers of derivatives of the gauge field j5|: 

1 - ±Re TrU pl = rj pl) £ Tr (F^) + a 2 [rft + r^R 2 + r^R 3 

" (1) 

+ 0(a 4 ) + total derivatives, 

R\ = Eju,!/ ^(D^F^D^F^), 
R2 = E^a TriD^D^), (2) 

Here the are coefficients in the operator product expansion of the plaquette. 
(Tree-level r^'s are tabulated for a variety of loop operators, like the plaquette, 
in 0.) Note that R\ communicates the lattice's violation of Lorentz 
invariance, and is therefore responsible for the bad behavior of the static 
potential, while R2 and R3 are Lorentz invariant. If we want to eliminate 
the 0(a 2 ) lattice artifacts then we need to form an improved action by 
adding other Wilson loops to the action, which will give canceling amounts 
of Ri, R 2 , R3. 

Only Ri contributes in Eq. (|l|) at tree-level, but quantum corrections 
bring in the other two operators. To remove all three, it might seem that we 
need to add three new Wilson loops to cancel these terms, but actually the 
coefficient of -R3 can be set to zero by a change of field variable in the path 
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Figure 1: Static-quark potential computed on 6 4 lattices with a ~ 0.4 fm 
using the f3 = 4.5 Wilson action and the improved action with f3 p \ = 6.8. 



integral, 



A^^A„ + a 2 a s f(a s ) £ D V F, 



(3) 



so only two new terms are needed [|10|, || . There are many possible choices 
for these, but we take the rectangle and "parallelogram" : 



rt 



17, 



pg 



(4) 



The improved action is [J 

= /? pl £§ReTr(l-[/ pl ) 
pi 

+ /3 rt £§ReTr(l-[/ rt ) 

rt 

+ /3 pg £±ReTr(l-£/ pg ), 



(5) 



pg 



with /3 p i given as an input, and /3 rt and /3 pg computed in tadpole-improved 
perturbation theory to cancel out the 0(a 2 ) terms in the derivative expansion 
of the action. At tree-level, the /3's are readily computed by combining 
expansions like Eq. (HD for each of the three loops. They are tadpole-improved 
by dividing each Wilson loop with L links by («o) L 0- One-loop corrections 
have also been computed fllTf , but must be adjusted to account for the tadpole 
improvement. We find: 



Prt 



20 Uq 



'1 + 0.4805 a & 



0.03325 a s . 

«0 



(6) 
(7) 



Following ||, we use the measured expectation value of the plaquette 
to determine both the value of the mean link ti and the QCD coupling 
constant a s , 



1/4 



u = (§ReTr([/ pl ) 

ln(|ReTr([/ pl ) 



3.06839 



(8) 
(9) 
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This result follows from the tree-level perturbative calculation of the 



plaquette |12| . The couplings /3 rt and (3 pg are determined self-consistently 
with uq and a s for a given f3 v \. As in NRQCD, there is no tuning of the 
couplings for the correction terms: tadpole-improved perturbation theory 
determines them in terms of the single bare coupling (3 V \. Using identities 
from || we find that our action is positive semidefinite at least for f3 p \ > 6.8, 
which is necessary if perturbation theory is to be reliable. 

The perturbative expansion for /3 rt already demonstrates the power 
of tadpole improvement. If one omits the tadpole factor Uq, the expression 
becomes /3 rt = — /3 p i(l + 2.0 a s )/20. Note that the coefficient of the one-loop- 
term a s has quadrupled. Tadpole improvement automatically supplies 75% 
of the one-loop contribution needed without improvement. Since a s m 0.3, 
the unimproved expansion is not particularly convergent. However, with 
tadpole improvement, the one-loop correction is only about 10-20% of (5 n . 

As indicated above, our improved action is not unique. Our techniques 
should work as well for actions with other forms for the correction terms. 
To verify this, we compare results obtained using our improved action above 
with those obtained from a tree-level improved action with a very different 
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s M [u] = /? P i£! ReTr (i-^ P i) 
pi 

+ /? trt £§ReTr(l-[/ trt ), (10) 



trt 

where U tT t is a twisted rectangle operator, 



u tTt = QQ • ( n ) 

At tree- level, with tadpole improvement, 

Note that tadpole improvement introduces four powers of uq here, rather 
than the two powers in /? rt above. This makes S tr t much more sensitive to 
tadpole improvement. 
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Action 


/3 P i 


(±ReTrf/ pl > 


aM c ° 


a 


dimensions 


improved 


6.8 


.46 


2.40 


.40 fm 


6 4 , 6 3 x 9 




7.1 


.506 


1.92 


.33 fm 


7 3 x 10 




7.4 


.56 


1.30 


.24 fm 


8 3 x 10, 9 3 x 12 


Strt 


4.1 


.454 


2.40 


.40 fm 


6 3 x 9 


Wilson 


4.5 


.34 


2.00 


.40 fm 


6 4 




5.7 


.55 


.80 


.17 fm 


12 3 x 24, 16 4 



Table 1: Parameters used in our simulations. The additional couplings for 
improved actions were determined using the values listed in the table for /3 pl 
and the plaquette. M? is the bare charm-quark mass used in the NRQCD 
quark action. 

3 Monte-Carlo Results 

We conducted simulations with improved actions for a range of large lattice 
spacings (Table [l]). The static-quark potential computed using our improved 
gluon action and the coarsest lattice is shown in Fig. |T]b. As in the Wilson 
case (Fig |l]a), the lattice spacing is about 0.4 fm. The dashed line in these 
plots is the standard infrared parameterization for the continuum potential, 
V(r) = Kr — 7r/12r + c, adjusted to fit the on-axis values of the potential. 
Off-axis points deviate from the fit by 40% for the Wilson theory, indicating 
a significant failure of rotation invariance due to fmite-a errors. By contrast, 
the deviations are only 2-4% for the improved theory — negligible for most 
low-energy applications. 

To assess the relative importance of tree- level improvement, tadpole 
improvement and one-loop corrections we computed the potential for several 
different actions, all with lattice spacings of about 0.4 fm. We focused on 
the deviation AV of V(a, a, a) from the continuum potential adjusted to 
fit on-axis values of the simulated potential. AV^ is a sensitive indicator of 
violations of rotational invariance. Our results are in Table |[ As expected, 
the correction term in the action is significantly underestimated without 
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Action 



AV{V3a)/KV3a 



unimproved (Wilson) 


.41 


(2) 


tree-level improved, no tadpole improvement 


.15 


(1) 


one-loop improved, no tadpole improvement 


.12 


(2) 


tree-level improved, with tadpole improvement 


.05 


(1) 


one-loop improved, with tadpole improvement 


.04 


(1) 


twisted-rectangle correction, with tadpole improvement 


.04 


(2) 



Table 2: Error in the static quark potential at V(a, a, a) for a variety of gluon 
actions. The lattice spacing in each case is a « 0.4 fm; K is the slope of the 
linear part of the static potential. 

tadpole improvement. The tadpole-improved action is very accurate both 
with and without one-loop corrections, suggesting that 0(a 2 a s ) corrections 
are comparable to those of 0(a A ). 

We have also included in Table |] results obtained using the twisted- 
rectangle action S trt ([[(]). We find that this action gives a potential that is 
essentially identical to that obtained with our other improved action. Note 
that tadpole improvement more than doubles the size of the correction term 
in Strt when a = 0.4 fm. The quality of the results obtained from S tv t is 
strong evidence in support of tadpole improvement. 

In computing the potentials, we assumed that the static quark propagator 
is simply the product of link operators U^(x) along the time axis. Our action 
is designed so that this is true at tree level through 0(a 2 ), however there are 
corrections of order a s a 2 . These arise because of the field transformation in 
Eq. (|3p. Our potential therefore has errors of 0(a s a 2 ), even though the action 
is accurate up to errors of 0(a 2 a 2 , a 4 ). Note that these particular errors in 
our potential do not break rotational invariance and so have no effect on 
the values of AV in Table |2|. A straightforward perturbative calculation is 
needed to remove the 0(a s a 2 ) errors from our potential. 
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Figure 2: S, P, and D states of charmonium computed on lattices with: 
a = 0.40 fm (improved action, /3 pl = 6.8); a = 0.33 fm (improved action, 
/3 p i = 7.1); a = 0.24 fm (improved action, /3 P \ = 7.4); and a = 0.17 fm 
(Wilson action, (3 = 5.7, from [6]). The dashed lines indicate the true masses. 

To further check on our improved theory, we examined the spin-averaged 
spectrum of the if) family of mesons using NRQCD for the c-quarks and 
our improved action at (5 V \ = 6.8, 7.1 and 7.4. Because we were examining 
only the spin-averaged spectrum, we omitted all spin-dependent corrections 
from the NRQCD action, but kept the corrections for 0(a, a 2 ) errors, and 
for spin-independent (9(v 2 /c 2 ) effects 0. The spectra, normalized to give 
the correct IP-IS splitting, are shown in Fig. ^together with experimental 
results (dashed lines) and simulation results obtained using the Wilson action 
for the gluons at smaller a's@. All agree to within 30-40 MeV. (The 
2S simulation masses are a little above the true mass, which is expected 
since the simulation does not include light-quark loops.) 

In Fig. H we show simulation results for the Coulomb-gauge radial 
wavefunctions of IS and IP charmonium as computed with the improved 
action on our coarsest and finest lattices. We also show wavefunctions 
computed from a continuum quark model tuned to reproduce our lattice 
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Figure 3: The radial wavefunctions for the IS and IP charmonium computed 
using improved actions and two different lattice spacings. Wavefunctions 
from a continuum quark model are also shown. Statistical errors are 
negligible for the IS* wavefunction. 
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Action 


/3 P i 


a(lP- IS) 


a 


(ip — is)/Vv' 


Wilson 


4.5 


1.410 (30) 


.41 fm 


1.38 (5) 




5.7 


.383 (10) 


.17 fm 


.87 (3) 


Improved 


6.8 


.922 ( 6) 


.40 fm 


.90 (1) 




7.1 


.756 ( 4) 


.33 fm 


.92 (1) 




7.4 


.558 ( 6) 


.24 fm 


.89 (2) 


Strt 


4.1 


.930 (30) 


.40 fm 


.92 (3) 



Table 3: Ratio of the charmonium IP — IS* splitting to the square root of 
the derivative of the static potential V(r) at 0.6 fm. Charmonium results for 
the Wilson action at f3 — 5.7 are from [6]. 



results. The wavefunctions from the different lattice spacings agree well 
everywhere except at r = 0, where agreement is not expected (because 
of renormalization effects). The wavefunctions show remarkable rotational 
symmetry considering that the rms radius of the 15 wavefunction, for 
example, is just 0.37 fm, or slightly less than one lattice spacing on the 
coarse lattice — accurate modeling of a hadron is possible with improved 
actions even when the hadron is only a few lattice spacings in extent. 

We also calculated the ratio of the IP-IS" charmonium splitting to the 
square root of the slope V of the static potential at r = 0.6 fm |14 



our 



results are in Table |3|. (We choose 0.6 fm because we have good simulation 
results for V(r) at that radius and, also, V is almost independent of r there.) 
This ratio is a dimensionless quantity that should become independent of 
lattice spacing as we approach the continuum limit. We see that, to within 
our errors, the improved action has reached the continuum limit at a lattice 
spacing of 0.4 fm, whereas the Wilson action has not. The twisted- rectangle 
action also gives excellent results at this lattice spacing. 

Since coarse lattices have far fewer sites and much less critical-slowing- 
down, the cost to produce a statistically independent configuration should 
be much less on a coarser lattice. To examine this issue, we compared our 
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results for V(r) with those in [[OJ which are for the potential computed 
using the Wilson action at (3 — 6 (32 4 lattice with a ~ 0.1 fm). We 
rescaled the coordinates and potential from this other study to put them 
in the same units as our (3 p \ = 6.8 results, and examined the potentials 
at comparable distances. The results from both simulations are listed in 
Table |j. In both cases the potential is obtained from the time dependence 
of loop-like correlation functions for times equal to or larger than some T min . 
Results should become independent of T min once it is sufficiently large; with 
our much smaller statistical errors, we can measure the small shift when T m ; n 
is doubled from .4 fm to .8 fm. Our results required 1.3 x 10 7 site updates, 
while the analysis on the fine lattice required 6.4 x 10 9 site updates. Since 
statistical errors (for T m j n = 0.4 fm) are about 20 times smaller for the coarse 
lattice, we estimate that comparable errors with the fine lattice would require 
197,000 times more site updates than we used on the coarse lattice. 

4 Conclusions 

We have found that, by using perturbatively-improved actions with tadpole- 
improved operators, we can accurately simulate quark and gluon dynamics 
on lattices with a spacing as large as 0.4 fm. Using either of two very different 
improved gluon actions, we obtained results for the confining potential and 
for charmonium that are independent of finite-lattice-spacing artifacts to 
within a few percent. 

It is striking that we can obtain accurate results on such coarse 
lattices using relatively simple actions. Following the ideas of the 
Wilson renormalization group, one expects, even at fixed, finite lattice 
spacing, to achieve arbitrarily accurate results by adding more and more 
operators to the discretized action, and by calculating their coefficients 
with increasing accuracy. We find that including just the leading 
corrections, with coefficients calculated in (tadpole-improved) tree-level 
perturbation theory, gives excellent results even on our coarsest lattices. 
Ultimately, nonperturbative calculations (for example, using Monte Carlo 
renormalization group methods) should be used to determine the coefficients, 
at least as a check of the perturbation theory. Such nonperturbative 
determinations of the couplings in the lattice action would be very interesting, 
although we expect such an effort to be much more difficult than the tadpole- 
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r/a c a = a c = .40 a = .10 

-^min -4 T m i n .8 T m [ n .3 

.962 .839 (11) 

1 .871 ( 0) .887 ( 1) 

1.361 1.316 (25) 

y/2 1.373 ( 1) 1.384 ( 2) 

1.667 1.581 (54) 

\/3 1.718 ( 2) 1.742 (10) 

1.924 1.842 (42) 

2 1.897 ( 2) 1.941 (10) 



Table 4: Comparison of the static-quark potential a V(r) as computed on a 
coarse lattice (improved action, /3 p i = 6.8, a = 0.40 fm), and on a fine lattice 
(Wilson action, /3 — 6, a — 0.10 fm). T min is the shortest time interval used 
in the correlation functions that determine V(r). All times and distances are 
in fm. 



improvement program outlined in this paper, and to produce only small 
changes to the coefficients determined here. In ||, the normalization of 
the strong coupling constant and the leading gluonic operator (F£ v ) was 
carefully studied. We found that tadpole-improved perturbation theory gave 
normalizations that were extremely close to the nonperturbatively obtained 
normalization. 

The coarseness of the lattice makes our simulations 10 4 -10 5 times faster 
than ones using unimproved actions; in fact, most of the a = 0.4 fm results 
in this paper were obtained using an IBM RS6000/250 desktop workstation, 
which is powered by a personal-computer CPU (66MHz PowerPC). We 
therefore believe that, using the methods we have described, the QCD 
hadron spectrum can be calculated to an accuracy of a few percent using 
computer resources that are already widely available. And by combining 
these techniques with forefront computing technology, we can begin to 
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tackle problems in nonperturbative QCD far more complex than previously 
imagined. 
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